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OPTIMIZATION APPLICATIONS 


The efficient design of aircraft structure involves a series of compromises 
among various engineering disciplines. These compromises are necessary to ensure 
the best overall design. To> effectively reconcile the various technical con- 
straints requires a number of design iterations, with the accompanying long elapsed 
time. Automated procedures can reduce the elapsed time, improve productivity and 
hold the promise of optimum designs which may be missed by batch processing. 

This presentation includes several examples of optimization applications 
including aeroelastic constraints. Particular attention is given to the success or 
failure of each example and the lessons learned. The specific applications are 
shown in Figure 1. The final two applications were made recently. 
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COPS ANALYTICAL MODEL 


Figure 2 illustrates the modeling of the stabilator in the Computerized 
Optimization Procedure for Stabilators (COPS); Reference 1 describes the procedure. 
The analytical model is a single-cell torque box idealized by eight discrete rigid 
chord streamwise sections with three mass points per section. Quasi -steady aero- 
dynamic forces act at user specified locations in each section. Nondimensional 
geometrical design parameters may be specified for taper ratio, thickness ratios 
at root and tip chords, aspect ratio, leading edge sweep angle, tip cut-off angle, 
pitch axis hinge line angle, pitch axis intersection with the mean aerodynamic 
chord (MAC), and spar locations. 
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Figure 2 
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COPS CONCEPTUAL FLOW DIAGRAM 

A greatly simplified flow diagram is shown in Figure 3. The procedure synthe- 
sizes, from the input data, a stabilator which satisfies all system constraints 
except those for flutter and divergence. A systematic perturbation of design vari- 
ables for 1) torsional stiffness, 2) balance weight, 3) pitch restraint and 4) roll 
restraint follows until the aeroelastic constraint is satisfied for minimum addi- 
tional weight. The procedure may be used in its basic sequential optimization 
scheme, where a new dynamic system is established after each iteration step. It 
may alternately be used in its simultaneous optimization mode, where each design 
variable is individually and exclusively evaluated from the same initial design 
point. The basic COPS program contains a realistic representation for every signi- 
ficant aspect of a believable stabilator flutter analysis and is fast enough, on 
the computer, to be used as an integral part of more encompassing aircraft systems 
optimization programs, as shown in the figure. 
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COPS EXAMPLE - EARLY CONCEPTUAL DESIGN FOR F-15 STABILATOR 

Using the semi-automatic "simultaneous" procedure, the COPS program was used to 
calculate the ratios of the change of flutter dynamic pressure to weight change 
(^/AW) for separate perturbations of stiffness at each of the elastic axis sta- 
tions, as shown in Figure 4a. Flutter was calculated for specified levels of AQ/AW 
* and compared, in Figure 4b, with optimization runs based on a torsional stiffness 
distribution proportional to the fourth power of the local chord and separately by 
balance weights at the tip leading edge. The balance weight of approximately 15 
lb is the minimum weight solution. Figure 4c shows the stiffness distributions 
for both the C^ and sensitivity approaches. 
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DETAIL DESIGN OF OPTIMUM F-15 STABILATOR 


Two separate configurations were considered for the final F-15 detail design, 
as shown in Figure 5a and discussed in Reference 2. The flutter model test results 
are summarized in Figure 5b. The 15-lb balance weight produces an overall increase 
in flutter speed with Mach number. The snag leading edge produces an overall 
increase in flutter speed, similar to that for the balance weight, at low speeds. 
However, the speed variation with Mach number is quite different, with the snag 
showing an initial sharper drop with increasing Mach number followed by a subse- 
quent sharper rise with further Mach number increase. Analyses indicated the 
favorable sharper rise to be associated primarily with the aft shift in stabilator 
aerodynamic center attributable to the area removed by the snag. The snag offered 
a significant weight savings over the balance weight with no effect on subsonic 
drag, aircraft stability or flying qualities. A small supersonic drag penalty was 
offset by the attendant weight reduction. 
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Figure 5 
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AEROELASTIC TAILORING STUDY CONFIGURATIONS 


Studies have been conducted on the use of the directional properties of compos- 
ite material to provide design improvements for fighter aircraft as discussed in 
Reference 3. The TSO (Aeroelastic Tailoring and Structural Optimization) computer 
program. Reference 4, which was developed by the Air Force Flight Pynamics Labora- 
tory (AFFDL), was used in these investigations. The configurations evaluated, 
shown in Figure 6, covered a wide spectrum of fighter aircraft aerodynamic 
surfaces, including 1) the F-15 composite wing, 2) a preliminary design horizontal 
tail, 3) a prototype aircraft movable outer panel, and 4) a conceptual wing for a 
future aircraft. The TSO program was validated with the F-15 composite wing which 
was designed to have the same distributed stiffness characteristics as the produc- 
tion metal wing. In spite of the structural approximations required by the TSO 

program, the predicted aeroelastic properties were surprisingly close to measured 
val ues. 
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AEROELASTIC TAILORING RESULTS 


Aeroelastic tailoring can play a significant role in the design of aircraft in 
various ways, as indicated in Figure 7. Specific detail is given for each configura- 
tion in References 3 and 5. 

As currently configured, the TSO computer program is appropriate for use pri- 
marily in preliminary design. The restrictive structural modeling requirements of 
TSO lead to converged results which are generally qualitative and which must be 
liberally interpreted when converting to a design that can be built. The experi- 
ence gained in the validation studies of the F-15 composite wing design, however, 
indicates that skillful use of the procedure can also yield good results in final 
detail design. 
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Figure 7 
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EQUIVALENT AFT SWEPT WING MODELS 


The optimized forward swept wing (FSW) was compared with three equivalent aft 
swept wings (ASW), shown in Figure 8, and evaluated for the same design constraints, as 
discussed in Reference 5. They are 1) an Equivalent Leading Edge sweep, where the 
ASW leading edge sweep angle is the negative of the FSW leading edge sweep angle, 
K 2) an Equivalent Elastic Axis sweep, and 3) a Flipped Wing. The wing geometry 

applies to all four wings. The wings were shifted longitudinally to give the same 
locations for the mean aerodynamic chords (MAC). The same aerodynamic and struc- 
tural models were used for all four wings. One of the apparent effects of the 
^ equivalent leading edge design is a structural bending axis that is about 20% 

shorter than the axis of the FSW. The bending axis for the flipped wing design, on 
the other hand, is about 12% longer. 
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COMPARISON OF FORWARD SWEPT WING WITH 
THREE EQUIVALENT AFT SWEPT WINGS 

Each of the ASWs was optimized by TSO and the results are shown in Figure 9. 
The FSW has the highest torque box skin weight and the Equivalent LE ASW, which is 
essentially a straight wing such as on the F-18, has the lowest weight. This 
weight advantage of the nearly straight wing is a direct result of the reduced 
structural axis length, which can be seen by comparing the bending moment normal to 
the elastic axis at the fuselage moldline. The air loading is also most favorable 
for the design of the fuselage carry-through structure on the FSW and the straight 
wing, as shown by considering both pitch and roll moments at the wing root. The 
ASWs are divergence free but have an active flutter constraint. The FSW has favor- 
able flutter properties, primarily because the frequency of the wing bending mode 
changes very little with increasing airspeed. Coupling with the torsion mode still 
occurs, but at a higher velocity than for the ASWs. 
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FASTOP APPLICATION TO NASTRAN BEAM-ROD 
DYNAMIC MODEL 


The FASTOP (Flutter and Strength Optimization Program) computer program (Refer- 
ence 6) which was developed by the AFFDL, has been applied to a beam-rod vibration 
and flutter idealization of a wing/store flutter model. The chosen configuration 
for this detail design application was a wind tunnel model with two stores on an 
outboard pylon and wing tip missile on, as described in Reference 7. The NASTRAN model 
is shown in Figure 10. The NASTRAN beam elements are based on GJ and El stiffness 
distributions, referred to an elastic axis, with similar distributions for the lead- 
ing and trailing edge control surfaces and the missile. There are rigid bars to 
connect the various components with the proper boundary conditions. Concentrated- 
elasticity members are used to represent integral springs, e.g. actuators, wing 
fold, missile/launcher/wing interfaces and wing/fuselage attachment. Structural 
optimization is not feasible because FASTOP does not calculate stresses in the beam 
elements. Steady air loads are not required because the starting point is an 
existing strength design. It was felt that the chances for success would be excel- 
lent for this simple straightforward model which has only 147 structural members. 
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FASTOP ANALYTICAL CONSIDERATIONS 


Many approximations and computational difficulties were encountered in convert- 
ing the NASTRAN model to FASTOP, as indicated in Figure 11. 
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• Vibration Calculated for Only 20 Normal Modes 

• Frequency Comparison Better Than Expected Considering Structural Compromises 

Unsteady Aerodynamics 
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Flutter Analysis 

• Non-Optimum Weight Factors Defined for Each Beam Element 
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FASTOP FLUTTER OPTIMIZATION RESULTS 


The results shown in Figure 12 look promising until one examines the redesign 
changes in the individual elements. The first 3 design cycles add increments of 
weight to the structural elements in proportion to their flutter velocity deriva- 
tives, 3Vf/8Wi, provided the derivatives are larger than an arbitrary minimim. 
This arbitrary minimum, which is not specified by the user, leads to an uneven span- 
wise distribution with peaks and valleys. It suffers from the lack of a built-in 
French Curve, which would smooth out the peaks and valleys to create a near-optimum 
design that could be built. The design cycles 4-10 continue the optimization by 
adjusting the weight distribution, while maintaining the desired flutter velocity. 
This weight adjustment reduces the increments along the wing torque box and builds 
up a large mass at the leading edge of the wing. The final design has only a large 
mass at the leading edge, near mid span, much like a forward mounted engine on a 
transport aircraft wing. 
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SENSITIVITY ANALYSIS DATA PREPARATION 

This approach to flutter optimization is based on a "sensitivity" technique 
similar to that first explored in conjunction with the development of the COPS 
program, which is described in Figure 4. The study was done in an extremely short 
elapsed time using existing flutter data sets for a beam-rod stabilator idealization 
based on NASTRAN and Doublet lattice. The only new data required are the GJ versus 
number of 45® plies per skin for several elastic axis (EA) stations, as shown in 
Figure 13a. With these data the change in GJ versus elastic axis station can be 
calculated for various weight increments, as shown in Figure 13b. 
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SENSITIVITY OPTIMIZATION RESULTS 


The steps in the sensitivity optimization study are given in Figure 14. Each 
step of the redesign is based on batch submittals of the NASTRAN flutter routines, 
followed by a conscious choice for the elements to be used in the subsequent step. 
After step 4 it is possible, by the use of engineering judgment, to specify a 
redesign distribution which satisfies the flutter requirement and is practical to 
build. These studies are state of the art in all respects and are quicker, cheaper 
and more accurate than possible with any currently available automatic optimization 
proced ure. 
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CONCLUSIONS AND RECOMMENDATIONS 


The following conclusions are appropriate in the area of flutter optimization. 

1) COPS, or a similar routine, is suitable for use in conceptual design for 
individual lifting surfaces when the geometry is undefined. 

2) TSO is suitable for use in preliminary design for individual lifting 
surfaces when the geometry is defined but there is still no well-defined 
structural model. Limitations in the structural model allow only limited 
use in detail design. 

3) FASTOP is based on a sound concept but is not general enough for use in 
detail design and is too difficult to use in preliminary design. 

4) NASTRAN based semi-automatic sensitivity techniques are the preferred 
approach for detail designs when the structural model is well defined but 
the flutter speed is deficient. 

Since NASTRAN is the accepted industry standard for structural analyses it 
seems appropriate to either 1) incorporate flutter optimization routines in NASTRAN 
or 2) ensure that any alternate program developments have a complete one-to-one rela- 
tionship to NASTRAN in all respects, and are designed to include generation of input 
data by graphics procedures. 


Conclusions 

• COPS Is Suitable for Conceptual Design 

• TSO Is Suitable for Preliminary Design 

• NASTRAN Sensitivity Technique Is Suitable for Detail Design 

Recommendations 

• Incorporate Flutter Optimization in NASTRAN 

• If Alternate Procedure^ Ensure Complete One-to-One 
Relationship to NASTRAN Including Graphics Generation of 
Bulk Data 


Figure 15 
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